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The continuum discretized coupled channels (CDCC) method is compared to the exact solution 
of the three-body Faddeev equations in momentum space. We present results for: i) elastic and 
breakup observables of d+ 12 C at Ed = 56 MeV, ii) elastic scattering of d+ 58 Ni at Ed = 80 MeV, and 
iii) elastic, breakup and transfer observables for 11 Be+p at Eu Be /A = 38.4 MeV. Our comparative 
studies show that, in the first two cases, the CDCC method is a good approximation to the full 
three-body Faddeev solution, but for the 11 Be exotic nucleus, depending on the observable or the 
kinematic regime, it may miss out some of the dynamic three-body effects that appear through the 
explicit coupling to the transfer channel. 

PACS numbers: 24.10.-i, 24.10.Eq, 25.55.Ci, 25.55.Hp, 25.60.Bx, 25.60Gc, 25.60.Je 



I. INTRODUCTION 

The strong coupling between elastic and breakup chan- 
nels in direct nuclear reactions involving deuterons led 
to the development of the continuum discretized coupled 
channels (CDCC) method where an effective three-body 
problem is solved approximately via the expansion of the 
full wave function in a selected set of continuum wave 
functions of a given pair subsystem Hamiltonian. Initial 
work by Johnson and Soper [1] showed that deuteron 
breakup was very important to understand reactions in- 
volving the deuteron. In that work, a two channel prob- 
lem was solved where the deuteron continuum was repre- 
sented by a single discrete s-state. Later developments by 
Rawitscher [2] and Austern [3] helped to introduce a more 
realistic representation of the continuum; further numer- 
ical implementations of the method proved its feasibil- 
ity [4] . Originally applied to reactions with the deuteron 
(e.g. [5]), it has since been extended to describe reactions 
with radioactive nuclear beams (e.g. [6-10]), namely to 
study elastic, transfer and breakup cross sections that re- 
sult from the collision of a halo nucleus with a proton or 
a stable heavier target such as 12 C or 208 Pb. 

In a recent paper [11] CDCC results obtained with two 
different basis sets, namely the basis set using the contin- 
uum of the projectile in the entrance channel, and the one 
using the continuum of the composite system in the final 
transfer channel, led to substantially different breakup 
cross sections for p( 11 Be, 10 Be)pn. These findings raise 
concern about the accuracy of the CDCC method, not 
only as a means to describe reaction dynamics, but also 
as an accurate tool to extract structure information on 
halo nuclei. 

An alternative approach to the solution of effective 
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three-body problems is the solution of the Faddeev equa- 
tions [12] for the wave function components or the equiva- 
lent Alt, Grassberger and Sandhas (AGS) equations [13] 
for the transition operators. The application of exact 
Faddeev/ AGS equations to the study of direct nuclear 
reactions has been shadowed in the past by the diffi- 
culty in dealing with the long range Coulomb force be- 
tween charged particles [14]. In recent calculations for 
the reaction c?+ 12 C [15] separable potentials are used 
and the Coulomb is taken into account only approxi- 
mately Given the progress achieved recently [16] for 
p-d elastic scattering and breakup we can now address 
the solution of effective three-body systems where two 
of the particles have charge. This was done first for d-a 
elastic scattering and breakup [17] and later for p- 11 Be 
elastic scattering and breakup [18] where n Be is a halo 
nucleus made up of a neutron and an inert 10 Be core. 
More recently the same system was used to study the 
convergence of the Faddeev/ AGS multiple scattering se- 
ries [19] at intermediate energy as a means to test the 
Glauber method. In all these works [16-19] we solve the 
AGS equations without resorting to a separable repre- 
sentation of the underlying interactions. Therefore the 
corresponding two-vector-variable integral equations are 
numerically solved without any approximations beyond 
the usual partial-wave decomposition and discretization 
of momentum meshes. 

Although CDCC was initially introduced as a practi- 
cal way of solving a complicated three-body scattering 
problem through a set of coupled Schrodinger-like equa- 
tions, later works [20, 21] tried to obtain a more formal 
justification of the method, by relating it to a trunca- 
tion of an orderly set of Faddeev equations. Further- 
more, it is argued that the CDCC solution approaches 
the exact (Faddeev) solution as the model space is in- 
creased. Although qualitative arguments are provided 
in those works to support the conclusions, they lack a 
numerical comparison between the CDCC and Faddeev 
methods in specific cases. The possibility of performing 
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this comparison, thanks to the recent developments in 
the numerical implementation of the AGS equations, is 
another motivation for the present work. 

Given these important new developments we propose 
here to benchmark, in a few test cases, CDCC results 
with exact solutions of the AGS equations. For this com- 
parison, we have selected the reactions d+ 12 C , d+ 58 Ni 
and n Be+p. The first two reactions correspond to classic 
cases for which there is data available [22-24] . The third 
reaction involves the scattering of a halo nucleus on a 
very light target, for which elastic, breakup and transfer 
have been measured before [8, 25, 26]. 

In Section II we present the AGS formalism, and in 
Section III we outline the different CDCC methodolo- 
gies we use. In Section IV we describe the details of the 
calculations and in Section V the results are presented. 
Conclusions are given in Section VI. 

II. THREE-BODY EQUATIONS 

This section provides the theoretical framework on 
which we base our calculations. Our treatment of the 
Coulomb interaction [16] relies on the screening and 
renormalization techniques proposed in Ref. [27] for two 
charged particle scattering and extended in Ref. [28] 
to three-particle scattering. The Coulomb potential is 
screened, standard scattering theory for short-range po- 
tentials is used, and the renormalization procedure is ap- 
plied to obtain the results for the unscreened limit. 

In the traditional odd-man-out notation of the three- 
body problem where pair (/?, 7) is denoted by a 
(a, (3, 7 = 1, 2, 3), the Coulomb potential w aR is screened 
around the separation r = R between two charged 
baryons (3 and 7. We choose w aR in configuration space 
as 

w aR {r)=w a (r)e-W> n , (1) 

where w a (r) — a e Zf)Z^/r represents the true Coulomb 
potential, with Zp (Z^) being the atomic number of par- 
ticle (3 (7), a e w 1/137 the fine structure constant, and 
n controlling the smoothness of the screening. We pre- 
fer to work with a sharper screening than the Yukawa 
screening (n = 1) of Ref. [14]. We want to ensure that 
the screened Coulomb potential w aR approximates well 
the true Coulomb potential w a for distances r < R and 
simultaneously vanishes rapidly for r > R, providing a 
comparatively fast convergence of the partial-wave ex- 
pansion. The screening functions for different n values 
are compared in Fig. 1, showing that the choice n = 4 
includes much more of the exact Coulomb potential at 
short distances than the Yukawa screening. In contrast, 
the sharp cutoff (n — > 00) yields an unpleasant oscillatory 
behavior in the momentum-space representation, leading 
to convergence problems. In Ref. [16] we found the values 
3 < n < 6 to provide a sufficiently smooth, but at the 
same time a sufficiently rapid screening around r = R; 
n = 4 is our choice in the present paper. 
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FIG. 1: Screening function wn(r)/w(r) as function of the dis- 
tance between the two charged particles r for characteristic 
values of the parameter n in Eq. (1): n = 1 (dashed-dotted 
curve) corresponds to Yukawa screening, n = 4 (solid curve) 
is the choice of this paper, and n — > 00 (dotted curve) corre- 
sponds to a sharp cutoff. 

We solve the AGS three-particle scattering equa- 
tions [13] in momentum space 

U$(Z) = 5 0a G o \Z) + J2~S^ R \Z)G (Z)U^(Z), 

(2a) 

U£\Z) = G \Z) +J2tW(Z)G (Z)u£\Z), (2b) 

(7 

where 5p a = 1 — 5p a , Gq{Z) is the free resolvent, and 
Ta R \z) the two-particle transition matrix derived from 
nuclear plus screened Coulomb potentials 

Ti R \Z) - {v a + w aR ) + {v a +w aR )G (Z)Ti R \Z), (3) 

embedded in three-body space. The operators U^{Z) 

and Uq R ^ (Z) are the three-particle transition operators 
for elastic/rearrangement and breakup scattering respec- 
tively; their dependence on the screening radius R is no- 
tationally indicated. On-shell matrix elements of the op- 
erators (2) between two- and three-body channel states 
0a(qi)^ai) and |</>o(p/q/)^0/) w hh discrete quantum 
numbers v ai , Jacobi momenta Pi and qi, energy E ail and 
Z = E a i + iO, do not have a R — > 00 limit. However, as 
demonstrated in Refs. [16, 28], the three-particle ampli- 
tudes can be decomposed into long-range and Coulomb- 
distorted short-range parts, where the quantities diverg- 
ing in that limit are of two-body nature, i.e., the on-shell 
transition matrix 

TaR (Z) = + W£ R G^ (Z)T^ (Z) , (4) 

Gi R \Z) = (Z-H -v a -w aR )- 1 , (5) 

derived from the screened Coulomb potential between 
spectator and the center of mass (cm.) of the bound 
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pair, the corresponding wave function, and the screened 
Coulomb wave function for the relative motion of two 
charged particles in the final breakup state. Those quan- 
tities, renormalized according to Refs. [16, 28], in the 
R — > oo limit converge to the two-body Coulomb scatter- 
ing amplitude (0 Q (q/)^ Q/ |T^-|</> Q (q;)^ Qi ) (in general, 
as a distribution) and to the corresponding Coulomb 
wave functions, respectively, thereby yielding the three- 
particle scattering amplitudes in the proper Coulomb 
limit 

(<t>0 (q/ ) vp t I U/3a | (f> a (q»)^ ai ) 

= 5 0a (<j) a (cif)v af \TZg-\<i> a (cii)is ai ) 

+ Jim {ZphQf){Mv>Pt\[u£\E ai + iO) 

- S Pa T^(E m + iO)]\M<k> at )KRMh 

(6a) 

_ l 

= i |^ o {^ 2 (P/)(0o(p/q/)^o / | (6b) 

x U^\E ai +i0)\<p a (q i )u ai )zj(q i )}. 

The rcnormalization factors Z aR (qi) and z R (pf) are di- 
verging phase factors given in Refs. [16, 27, 28]. 

Z aR (q) = e- 2iS «*^, (7a) 

where 5 a R{q), though independent of the relative angular 
momentum I in the infinite R limit, is realized by 

S aR (q) = of(q) - (7b) 

with the diverging screened Coulomb phase shift rj" R (q) 
corresponding to standard boundary conditions and the 
proper Coulomb one af(q) referring to the logarithmi- 
cally distorted Coulomb boundary conditions in channel 
a with orbital angular momentum / between particle-pair 
a. For the screened Coulomb potential of Eq. (1) the in- 
finite R limit of 5 aR (q) is known analytically 

5 aR (q)=K a (q)[ln(2qR)-C/n], (7c) 

K a (q) = a e Z a (Zp + Zj)M a /q being the Coulomb pa- 
rameter, M a the reduced mass, and C ~ 0.5772156649 
the Euler number. Likewise 

z R (p)=e- 2iS *to, (8a) 

where 

5 R (p) = K(p)[\n{2pR)-C/n) (8b) 

with n(p) = a e Zp Z^Ha/p where [3 and 7 denote the 
two charged particles and [i a their respective reduced 
mass. The R — > 00 limit in Eqs. (6) has to be calcu- 
lated numerically, but due to the short-range nature of 
the corresponding operators it is reached with sufficient 
accuracy at rather modest R if the form of the screened 



Coulomb potential has been chosen successfully as dis- 
cussed above. More details on the practical implementa- 
tion of the screening and renormalization approach are 
given in Ref. [16]. 

The three-body results are obtained from the solu- 
tion of the AGS equations (2) for nuclear plus screened 
Coulomb interaction together with the renormalization 
procedure (6). The equations are solved using partial 
wave decomposition and retaining as many channels as 
needed for convergence. Our numerical technique for 
solving AGS equations with non-separable potentials is 
explained in more detail in Refs. [29, 30] in the context 
of nucleon-deuteron scattering. 



III. CDCC FORMALISM 

The CDCC method [3, 4] was introduced as an approx- 
imate solution to the three-particle Schrodinger equa- 
tion. Its main objective is to provide a reliable yet prac- 
tical way of describing reactions involving three-body 
breakup. 

Let us consider specifically the breakup reaction p+t — > 
c+x+t. In CDCC, the wavefunction is expanded in terms 
of only one Jacobi coordinate set (r,R), 

POO 

*^(r, R) = £ Mr)^p(R) + / dkMrWk(R), 
P J o 

(9) 

where <j> p (r) are eigenfunctions of the projectile, p(k) be- 
ing a general subscript for projectile bound (continuum) 
states, and tjj^(R) the spectator wavefunction for the 
motion of the projectile relative to the target. The Ja- 
cobi coordinate r describes the c+x relative motion while 
R the p + t relative motion. In CDCC the Schrodinger 
equation expressed in this Jacobi set reads 

{H 3b -E)* K {r,R)=0. (10) 

where the three-body Hamiltonian is separated into the 
internal Hamiltonian of the projectile and the relative 
motion between the projectile and the target: H 3 b = 
H^t + T R + U xt + U ct , where H int = T r + V xc (r). The 
projectile is modelled by a real potential which produces 
its initial bound state, whereas the fragment-target in- 
teractions should contain absorption from channels not 
included explicitly in the model (optical potentials). In 
principle, this is important for the validity of the CDCC 
method [31], since it reduces the coupling to the other 
three-body channels best described by other Jacobi sets 
taken explicitly into account by the Faddeev method. 

For practical reasons, the integral over projectile scat- 
tering states (9) is discretized and truncated at a maxi- 
mum energy. There are several methods of discretization 
but here we use the average method, where the c + x 
scattering radial functions Uk (r) are averaged over k to be 
made square integrable [6, 11]. Thus, the radial functions 
for the continuum bins in the average method, u p (r), are 
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a superposition of the projectile scattering eigenstates 
I 2 r k P 

M r ) = J^j-J k g P (k)u k (r)dk, (11) 

with weight function g p (k). The normalisation constant 
is defined by N p = j k P \g P (k)\ 2 dk. 

After a few steps of algebra and using the eigenvalue 
equation for the projectile H int 4> P = s P ^> P , the standard 
CDCC equation becomes 

[T R +V PP (R)-E P ]^ (R) = -J2 V pp ,(R)iP*(R), (12) 

where E p — E am , — e p , and the coupling potentials con- 
tain both nuclear and Coulomb parts V pp i (R) — (<f) P \U x t+ 
U c t\4>p')- This equation is indeed a coupled channel equa- 
tion, coupling the projectile ground state to its contin- 
uum states via Vo p , but also coupling projectile states 
within the continuum, called the continuum-continuum 
couplings. The solution of the coupled equations provide 
the wavefunctions ip^(R). The scattering observables 
(associated to the elastic and breakup channels) are ex- 
tracted from the asymptotic behaviour of ipjf(R). Nu- 
merical solutions of Eq. (12) involve partial wave expan- 
sions of (j> p {r), ip^{R) and a multipole decomposition of 
Uxt + U ct . 

In its standard form, CDCC models the breakup of a 
projectile as inelastic excitation (CDCC-BU). However, 
the CDCC wavefunction of Eq. (9) can also be used in 
the exit channel of a transfer reaction. In that case, 
the breakup process is understood as a transfer of the 
fragment to the continuum of the final composite system 
(CDCC-TR*). In the CDCC-TR* scheme, the scattering 
observables can be obtained inserting the CDCC wave- 
function in the prior form of the transition amplitude: 

T pr ior Vxt + tfct-EWte), (13) 

where \ P is a distorted wave, generated by the poten- 
tial U p t(R)- In this work, we chose U p t{R) = (<f>o\U x t + 
Uct\<Po) as the single folding of the core-target and the 
fragment-target interactions over the projectile's ground 
state (the Watanabe potential). We notice that, if \fj ^ 
is the exact solution of the three-body Hamiltonian, the 
transition amplitude (13) is exact and does not depend on 
the choice of the auxiliary potential U p t- In practice, this 
wavefunction is replaced by an approximate one which, in 
the TR* approach, corresponds to the CDCC expansion 
in the final channel. 

The transfer to the continuum approach seems to 
provide a good description of the data in some cases 
[10, 32, 33]. In Ref. [11], a detailed discussion of these 
two mechanisms is presented along with a comparative 
study. In principle, and as long as the CDCC model 
space is sufficiently large, one would expect that differ- 
ent choices of the Jacobi coordinate produce the same 
results. However, in [11] it is shown that this equivalence 
docs not hold for a number of reaction observables. 



IV. DETAILS OF THE CALCULATIONS 

As mentioned above we study the scattering of 
deuterons on 12 C at E d = 56 MeV and 58 Ni at E d = 
80 MeV, as well as n Be on protons at E^^/A = 
38.4 MeV. All reactions are considered as effective three- 
body problems, namely p + n + 12 C, p + n + 58 Ni, and 
10 Be + n + p, where 12 C, 58 Ni and 10 Be are taken as 
inert cores. Therefore in the present section we de- 
fine the interactions between all pairs together with the 
model space used in solving the AGS and the CDCC 
equations. For simplicity all interactions are taken spin- 
independent as often done in many CDCC calculations 
(e.g. Ref. [3]) and therefore all particles are considered 
spinless bosons. Nevertheless, spin-independent interac- 
tions, especially for the n-p system, are only semirealistic, 
but are sufficient for a benchmark comparison. This is 
not a limitation of the CDCC method nor of the AGS 
method as we have already demonstrated in nuclear re- 
action calculations [17-19] with realistic spin-dependent 
potentials. 

A. d+ 12 C at E d = 56 MeV 

The interactions between neutron- 12 C and proton- 12 C 
are optical potentials that fit the elastic scattering at half 
the incident laboratory energy [34]; the parameters are 
taken from the global fit by Watson et al. [35]. The 
neutron-proton bound and continuum states are mod- 
elled with a simple Gaussian interaction fitted to the 
deuteron binding energy 

V(r) - -V e-( r / r °) 2 , (14) 

where V = 72.15 MeV and r = 1.484 fm. The same 
interaction is used in all three test cases and corresponds 
to the choice in [5, 37]. 

The model space needed for converged solutions of the 
AGS equations contains partial waves I < 3 in the n-p 
relative motion, I < 12 in the n- 12 C channel, and I < 24 
for elastic scattering (I < 32 for breakup) in the p- 12 C 
channel. This last channel is more demanding due to the 
presence of the Coulomb force. Total angular momen- 
tum up to J = 30 (J = 60 for breakup) is included. The 
Coulomb potential is screened with a radius of R = 10 fm 
(R = 18 fm for breakup) and smoothness n — 4 (see 
Eq. (1)). The exception is the breakup kinematical situ- 
ations characterized by small momentum transfer in the 
p- 12 C subsystem which are sensitive to the Coulomb in- 
teraction at larger distances and therefore need larger 
screening radius and a special treatment as described in 
the Appendix. Note that, as optical potentials are used 
for n- 12 C and p- 12 C, there is no transfer to p- 13 C and 
n- 13 N channels. 

The corresponding CDCC calculations include n-p par- 
tial waves I < 8 and bins up to £ max = 46 MeV (with 
15 (10) energy bins for the even (odd) partial waves, 
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evenly spaced in linear momentum) integrated up to 
^bin = 80 fm; for the total angular momentum we in- 
clude J < 60. The coupling potentials were expanded in 
multipoles (Q) up to Q maK = 6. We note that this rela- 
tively large model space is required to achieve sufficient 
accuracy for the breakup observables. If only elastic scat- 
tering is required, / < 2 gives almost a converged result. 
The CDCC equations are solved up to i? max = 100 fm. 

B. d+ 58 Ni at E d = 80 MeV 

Similarly to the Carbon test case, we study the scatter- 
ing of deuterons from a Ni target using neutron- 58 Ni and 
proton- 58 Ni optical potentials from the global parameter- 
ization of Becchetti and Grenless [38] , evaluated at half 
the incident laboratory energy. The n-p interaction is 
given by Eq. (14). 

In this case, the model space needed for converged solu- 
tions of the AGS equations for elastic scattering contains 
partial waves I < 3 in the n-p relative motion, I < 14 
in the n- 58 Ni channel, and I < 32 in the p- 58 Ni channel. 
Again this last channel is more demanding due to the 
presence of the Coulomb force. Total angular momen- 
tum up to J = 60 is included. The Coulomb potential 
is screened with a radius of R = 10 fm and n = 4 (see 
Eq. (1)). 

For the CDCC calculations we include n-p partial 
waves I < 2. For I = 0, 1 the continuum was truncated at 
E max = 30 MeV, and divided into 12 bins evenly spaced 
in the linear momentum, whereas for I = 2 we include 
excitation energies up to E max = 50 MeV, and use 20 
bins; Q < 2 multipoles are retained in the expansion of 
the coupling potentials. The coupled equations are inte- 
grated up to R max = 80 fm with total angular momentum 
up to J = 100. 

C. "Be+p 

To describe the scattering of 11 Be from a proton tar- 
get, we need a binding potential for the n- 10 Be pair, as 
well as fragment-target optical potentials. The n- 10 Be 
interaction takes the standard Woods-Saxon form 

V(r) = -Vof(r,R ,a ) (15) 

with 

/(r,i2,a)=(l + eC-*>/°)~ 1 ) (16) 

where Ri = riA^ and A the mass number of 10 Be. 
The geometry of the interaction is fixed, with a radius 
ro = 1.39 fm and a diffuscncss clq = 0.52 fm. The depth 
of the interaction is L-dependent and the corresponding 
values of V for each partial wave are given in Table I, 
together with the energies for the corresponding bound 
states and resonance. In the three-body calculation the 



lowest (Pauli forbiden) bound state |6o) in L = is moved 
to a large positive energy T, replacing the potential V by 
V' = V + \bo)T(bo\. In the r — > oo limit this is equivalent 
to projecting |6 ) out as demonstrated in Ref. [39]. In 
practical calculations we found that r « 2 GeV is suf- 
ficiently large in order to obtain T- independent results. 
Thus, the state with eg = -0.503 MeV is left as the 
ground state of 11 Be. The unphysical deep s-state is also 
left out of the CDCC calculations. 



TABLE I: Parameters of the n- Be interaction in different 
partial waves and resulting energies of bound states and res- 
onance. 



L 


Vo (MeV) 


e L (MeV) 


e* L (MeV) 





51.639 


-30.28 


-0.503 


1 


26.264 


-0.183 




2 


51.639 


1.317 -i 0.188/2 




> 3 


51.639 







As for the projectile-target interactions, given such a 
light proton target, one of the potentials is simply the n-p 
interaction that binds the deuteron and that is used in 
IV A and IV B (no absorption). The p- 10 Be is obtained 
from a direct fit to elastic data [25]. We use the standard 
optical potential form 

V(r) = -Vo f(r, Ro, o ) - i W v f{r, R v ,a v ), (17) 

where f(r,R,a) is given by Eq. (16), plus the Coulomb 
interaction of a uniform charge sphere with radius R c = 
r c Az . A good fit to the p- 10 Be data at E-^^/A = 
39.1 MeV and up to 6 C m = 70° (see Fig. 7) is ob- 
tained with the parameter set: V = 51.2 MeV, W v = 
19.5 MeV, r c = r = r v = 1.114 fm, a = 0.57 fm and 
a v = 0.50 fm, as used in [11]. These parameters are 
slightly different from the ones proposed by Watson et 
al. [35]. Since the energy is sufficiently close to 11 Be-p 
scattering at E^^/A = 38.4 MeV we expect the fit to be 
appropriate. 

The Faddeev model space contains partial waves I < 4 
in the n-p relative motion, / < 5 in the n- 10 Be channel, 
and I < 22 in the p- 10 Be channel. Again, this last channel 
is more demanding due to the presence of the Coulomb 
force. Total angular momentum up to J = 20 (J = 40 for 
breakup) is included. The Coulomb potential is screened 
with a radius of R = 10 fm and smoothness n = 4. An 
additional difficulty is the presence of the sharp d-w&ve 
n- 10 Be resonance which is treated using the subtraction 
technique as in Ref. [17]. 

CDCC calculations are performed using both the n Be 
breakup states (CDCC-BU), where the reaction mecha- 
nism is inelastic excitation of the projectile into its con- 
tinuum, and the deuteron breakup states (CDCC-TR*), 
where the reaction mechanism involves transfer to the 
continuum of the deuteron in the d+ 10 Be transfer chan- 
nel [11]. In the course of the calculations, it became ap- 
parent that the model space used in [11] was not enough 
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to achieve full convergence of the CDCC-BU calcula- 
tions. In this work, we have increased the number of 
partial waves for the n- 10 Be relative motion, as well as 
the number of multipoles for the coupling potentials, up 
to l max — 8, Q m ax = 8, respectively. Even with this 
large number of partial waves, the results were not com- 
pletely converged. Inclusion of higher partial waves led 
to numerical instabilities in the calculations, and hence 
the results presented here correspond to l max = 8. Con- 
tinuum bins were calculated up to E max = 34 MeV for 
I < 6 and E max = 32 MeV for 1 = 7,8. 

For CDCC-TR*, the model space is also augmented 
with respect to the calculations performed in [11]. The 
number of partial waves for the pn relative motion is 
increased to I < 8, and bins are considered up to E max = 
35 MeV and i?bin = 60 fm. We notice also that in the 
present CDCC-TR* calculations the n- 10 Be interaction 
in the final channel (d+ 10 Be) is real, while in [11] this 
interaction was complex. Multipoles Q < 4 are included 
for the CDCC coupling potentials. An extended non- 
locality range of 14 fm was used for the transfer couplings. 

For both the CDCC-BU and CDCC-TR* calculations 
the total angular momentum is J < 35, and the coupled 
equations are integrated up to i? max = 60 fm. 

V. RESULTS 

Our three test cases are chosen to span a variety of 
situations. The d+ 12 C and d+ 58 Ni reactions at inter- 
mediate deuteron energies contain important effects that 
could not be well accounted for by simple prior form 
DWBA calculations [3, 37]. These reactions have been 
measured before and detailed breakup data are available 
for the former [37]. Finally we include a reaction involv- 
ing a loosely bound halo nucleus, the study of 11 Be+p, 
where three-body breakup has a decisive contribution 
and for which there have been several experimental stud- 
ies (e.g. [25]). Previous studies [40] have suggested the 
relevance of the interplay betweeen the breakup and the 
(p, d) transfer channel. The possibility of including both 
channels within the Faddeev formalism constitutes a fur- 
ther motivation of the present analysis. 



A. d+ 12 C at E d = 56 MeV 

Our results for <i+ 12 C elastic scattering at Ed = 
56 MeV are presented in Fig. 2, and compared to data 
around the same energy. The differential cross section 
da/d^l is divided by the corresponding Rutherford cross 
section dan/dSl. First we point out that the CDCC cal- 
culation (dash-dotted line) reproduces the exact three- 
body results (solid line) up to very large scattering an- 
gles. Secondly, there is agreement with the data at for- 
ward angles (9qm < 60°) but this agreement deteriorates 
for backward angles where mechanisms other than three- 
body breakup may start to play a role. 



100 b j 
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FIG. 2: (Color online) Elastic cross section for deuterons on 
12 C at E c i — 56 MeV: the solid line corresponds to exact 
three-body results and the dash-dotted line to CDCC. The 
result of a single channel cluster folding (dotted) is also shown. 
The experimental data (diamonds) are from Ref. [22]. 
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FIG. 3: (Color online) Semi-inclusive differential cross section 
versus proton scattering angle for the breakup of deuterons 
on 12 C at Ed = 56 MeV: the solid line corresponds to the 
exact three-body results and the dash-dotted line to CDCC. 
The dotted line corresponds to the three-body exact result in 
the absence of the Coulomb force. The experimental data are 
from Ref. [37]. 



In order to show the well known influence of the 
deuteron continuum on the elastic channel we have also 
included in Fig. 2 the CDCC calculation without any 
coupling to the continuum (dotted line). This corre- 
sponds to a one-channel calculation with the deuteron- 
target potential given by the single-folding expression: 
Vqo{R) — {4>d\V p t + V n t\4>d)- It can be seen that this 
calculation largely overestimates the data in the angu- 
lar region where the cross section is large, evidencing the 
importance of the deuteron breakup channel in the dy- 
namics of the reaction. 

Next we consider the breakup observables. In [37], 
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measurements were taken by fixing the neutron detector 
at 9 n = 15 deg. We present both the proton angular dis- 
tribution (after integration over energy) and the energy 
distributions for specific proton angles and compare to 
the data. The agreement between CDCC (dash-dotted 
line) and the exact three-body results (solid line) is seen 
over all proton angles as shown in Fig. 3. Also shown 
in Fig. 3 are the three-body results obtained without 
the Coulomb interaction (dotted line), where it becomes 
clear that Coulomb cannot be neglected for a wide an- 
gular range at forward scattering angles. We note that 
the agreement of present CDCC calculations with the 
data has been improved compared to the CDCC studies 
presented in [5] , probably due to the larger model space 
included in our work. 

In addition a note of caution needs to be added vis-a- 
vis the convergence of the three-body results with screen- 
ing radius R for —30° < 9 p < 10° where we face the 
most demanding phase space constraints. Since this re- 
gion corresponds to small momentum transfer in p- 12 C 
subsystem, the convergence with screening radius is slow 
and not uniform, forcing us to use a more sophisticated 
treatment for the breakup amplitude as described in the 
Appendix. For this reason the accuracy of our calcula- 
tion is 10% to 15% for -30° < 9 p < 10° and better than 
5% for all other values of 9 p . The most sensitive region 
is the maximum of the cross section. 

The proton energy distributions are shown in Fig. 4 
for protons being scattered to the same side as the neu- 
tron (8 p > 0), and in Fig. 5 for protons coming out at 
opposite angles from the neutron (9 P < 0). The overall 
agreement with the data is remarkable. Missing cross 
section is visible for the large positive scattering angles 
starting with 8 p — 25°. As for the angular distribution, 
our CDCC results show an improvement compared to 
the analysis of [5] , which we attribute to the inclusion of 
I = 4 in the relative motion of n-p subsystem in our model 
space. Most important for this work is the realization 
that CDCC simulates the three-body effects contained 
in the solution of the exact three-body problem, even at 
this level of detailed observables. Again, for the reasons 
mentioned above, the three-body results shown in Fig. 5 
for 9 p — —15°, —20° and —25° may change slightly with 
the chosen screening radius while at all other angles we 
have fully converged results. 



B. d+ 58 Ni at E d = 80 MeV 

We have reproduced the elastic scattering results pre- 
sented in [3] , and compare in Fig. 6 our theoretical predic- 
tions to two different sets of data around the same energy. 
As before, the solid line corresponds to the three-body re- 
sults from the solution of the AGS equations, whereas the 
dash-dotted line is calculated using CDCC. Both calcu- 
lations agree perfectly up to scattering angles 9 d = 80°, 
providing a good description of the data. We also show 
the result with the cluster folding potential (dotted line) 



which demonstrates the importance of deuteron breakup 
in the reaction mechanism, as was the case for the elastic 
scattering of d+ 12 C. The slight disagreement of the full 
calculations and the data could be due to the ambiguities 
of the optical potentials, or due to the influence of other 
channels not included explicitly in our calculations (such 
as transfer or target excitation). 

C. ^Be+p at £ Lab /A=38.4 MeV 

Our last test case involves the breakup of the loosely 
bound 11 Be on a very light target, the proton. Previous 
works have found difficulties in describing this process 
[41] and, furthermore, this reaction raised the red flag 
when comparing two different CDCC calculations which 
should produce the same results [11]. We revisit the topic 
in the hope that the exact three-body calculations can 
help shed light on the issue. 

In Fig. 7 we show the results of our calculations for 
11 Be+p elastic scattering together with the correspond- 
ing data. For comparison we include p- 10 Be elastic data 
at E-Lab/A = 39.1 MeV and the corresponding theoretical 
fit (dashed curve) obtained with the p- 10 Be optical poten- 
tial given in Sec. IV C. Two important features immedi- 
ately arise: i) the agreement between the CDCC and the 
exact three-body results, and ii) the mismatch with the 
data. In this work the first point is of more relevance than 
the second, demonstrating that the CDCC takes well into 
account the three-body effects fully present in the AGS 
approach. However, point ii) suggests that in this reac- 
tion, degrees of freedom beyond three-body breakup are 
being excited [41]. 

One main difference between this and the previous two 
examples is the explicit inclusion of the transfer channel 
in the three-body calculations. In other words, there is 
no absorption in n-p while the corresponding interactions 
in the previous two test cases n- 12 C and n- 58 Ni included 
absorption. In u Be+p, the neutron transfer channel is 
very important. We show in Fig. 8 the three-body pre- 
dictions for the transfer 11 Be(p, rf) 10 Be (thick solid line) 
together with the data for E p = 35.5 MeV [26]. The 
three-body calculation predicts the transfer cross section 
~ 20% above the data. If a simple proportionality of 
the transfer cross section to the square of the n- 10 Be 
I = single particle wavefunction were to be assumed, 
the three-body results would suggest ground state spec- 
troscopic factors consistent with previous works [40]. 

We also show in this figure the prediction of the CDCC- 
TR* (dashed line), obtained with Eq. (13). Due to 
the impossibility of including partial-wave dependent in- 
teractions in the evaluation of the coupling potentials, 
this calculation was performed assuming n- 10 Be poten- 
tial (15) with Vo — 51.639 MeV in all partial-waves that 
reproduces the ground state of 11 Be. For a meaningful 
comparison with the exact result, we include also in this 
figure a Faddeev calculation performed with the same 
n- 10 Be interaction (thin solid line). Notice that in this 
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FIG. 4: (Color online) Exclusive differential cross section versus proton energy for the breakup of deuterons on 12 C at Ed = 
56 MeV, 8 n — 15° and 9 P > 0: the solid line corresponds to exact three-body results and the dash-dotted line to CDCC. The 
experimental data are from Ref. [37]. 
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FIG. 5: (Color online) Exclusive differential cross section versus proton energy for the breakup of deuterons on 12 C at Ed = 
56 MeV, for 6 n = 15° and 9 P < 0: the solid line corresponds to exact three-body results and the dash-dotted line to CDCC. 
The experimental data are from Ref. [37]. 



case there is no bound excited state in 11 Be. We see 
that the cross section for this CDCC-TR* calculation 
is about 15% smaller than the AGS, and would hold a 



spectroscopy factor closer to unity. The difference be- 
tween AGS and CDCC could be due to the fact that the 
CDCC wavefunction is not a good reproduction of the 
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FIG. 6: (Color online) Elastic cross section for deuterons on 
58 Ni at Ed — 80 MeV : the solid line corresponds to ex- 
act three-body results and the dash-dotted line to CDCC re- 
sults. The experimental data at 80 MeV (diamonds) are from 
Ref. [23], and those at 79 MeV (triangles) are from Ref. [24]. 
The dotted line is the CDCC calculation suppressing the cou- 
pling to the deuteron continuum (see text). 
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FIG. 7: (Color online) 1 H( 11 Be, n Be)p elastic cross section 
at i?Lab/j4 = 38.4 MeV. The solid line corresponds to exact 
three-body results while the dash-dotted line to CDCC. The 
dashed line corresponds to an optical potential fit to the cor- 
responding 10 Be-p data of Ref. [25] shown by the triangles. 
The diamonds correspond to n Be-p elastic data of Ref. [25]. 



exact three-body wavefunction in the surface region, or 
that the choice of the optical potentials appearing in the 
remnant term of Eq. (13) is inadequate for this purpose, 
which could be connected to the poor description of the 
11 Be elastic data. 

Finally in Figs. 9 and 10 we show the semi-inclusive 
differential cross section for the breakup 11 Be + p — » 
10 Be + p + n where 10 Be is the detected particle. We 
present both, the energy distribution (Fig. 9) and the an- 
gular distribution (Fig. 10). For the energy distribution, 




FIG. 8: (Color online) Transfer reaction ^("Be, 10 Be)d cross 
section at E^b/A = 38.4 MeV. The thick solid line cor- 
responds to exact three-body result, while the dotted line 
corresponds to the same calculation multiplied by 0.7. The 
thin solid line is the exact calculation with a partial-wave in- 
dependent n- 10 Be interaction. The latter is to be compared 
with the CDCC-TR* calculation (dashed line), as explained 
in the text. The experimental data are from Ref. [40] at 
E p = 35.3 MeV. 



two CDCC-BU calculations are shown, one with I < 8 
(dash-dotted line) and one with I < 6 (thin solid line) for 
the n- 10 Be motion. The significant difference between 
these two calculations suggests that the CDCC-BU cal- 
culation is not converged with respect to the number of 
ri- 10 Be partial waves. The calculation with I < 8 repro- 
duces reasonably well the shape of the energy distribution 
predicted by the AGS calculation, but it underestimates 
this cross section at the peak by about 20%. This un- 
derestimation could be due to the contribution of higher 
n- 10 Be partial waves or due to some breakdown of CDCC. 
The angular distribution is also in good agreement with 
the exact result for the whole angular range, except at 
very small angles. 

The CDCC-TR* calculation also reproduces reason- 
ably well the energy distribution. For the higher 10 Be 
energies, this calculation is however well above the AGS 
result. Moreover, some underestimation of the cross sec- 
tion is observed at the maximum of the distribution, as 
well as at low 10 Be energies. We notice that these ener- 
gies are associated with configurations in which the pn 
system is in a very high excited state, and these states 
are difficult to include in the CDCC-TR* calculation. 
The peak observed in the energy distribution of Fig. 9 
corresponds to n-p quasi-free scattering and it is nat- 
ural that it is best reproduced by CDCC-TR*. Note 
that, as for the CDCC-TR* predictions for the transfer 
to the ground state, in the CDCC-TR* breakup, the fixed 
ri- 10 Be partial-wave independent interaction was used to 
generate the deuteron continuum. However in this case 
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FIG. 9: (Color online) Semi-inclusive differential cross section 
for the reaction ^("Be, 10 Be)pn, at £ Lab /4 = 38.4 MeV, 
versus 10 Be center of mass energy. The thick solid line corre- 
sponds to exact three-body results, the dashed line to CDCC- 
TR*, the dash-dotted line and the thin solid line to CDCC-BU 
with Imax = 8 and l ma x = 6, respectively. 
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FIG. 10: (Color online) Semi-inclusive differential cross 
section for the reaction 1 H( 11 Be, 10 Be)pn, at E^^/A — 
38.4 MeV, versus 10 Be angular distribution after energy in- 
tegration. The solid line corresponds to exact three-body re- 
sults, the dashed line to CDCC-TR* and the dash-dotted line 
and the thin solid line to CDCC-BU, with lmax = 8 and 
lmax = 6, respectively. 



the results are not very sensitive to this potential choice. 

The CDCC-TR* angular distribution (dashed line in 
Fig. 10) reproduces well the Faddeev calculation at small 
angles, but underestimates the cross section for angles be- 
yond 50°. We notice again that small angles are mainly 
associated with the np quasi-free scattering region which 
is better described in a pn basis, as it is done in the 
TR* approach. Conversely, these configurations are dif- 
ficult to describe in the n- 10 Be basis which explains the 
low convergence rate of the BU calculation at small 10 Be 
scattering angles. 

We finish this section by noting that, even if the CDCC 



wavefunction is not appropriate to describe breakup in all 
regions of phase-space, this wavefunction can be accurate 
in a limited domain. For example, in the CDCC-BU 
calculations presented in this section, the CDCC solution 
can be a good approximation to the exact three-body 
wavefunction in the region of space that corresponds to 
small neutron- 10 Be separations. 



VI. CONCLUSIONS 

A comparative study of reaction observables calculated 
within the three-body AGS framework and the approx- 
imate CDCC equations is presented. The AGS results 
shown here involve heavier nuclei where optical potentials 
together with the full treatment of the Coulomb interac- 
tion are used to describe direct nuclear reactions driven 
by deuterons and halo nuclei. We perform calculations 
for the scattering of deuterons on 12 C at Ed = 56 MeV 
and 58 Ni at E d — 80 MeV, as well as 11 Be on protons at 
Ehab — 38.4 MeV /A, and calculate elastic, breakup and 
transfer observables. 

The results indicate that, for reactions involving the 
elastic scattering and breakup of deuterons on carbon 
and Ni targets, CDCC is in agreement with the full 
three-body results. Our calculations also reveal that 
Coulomb effects in the breakup of deuteron by 12 C are 
not negligible for proton forward-angle breakup. Indeed, 
the method of screening and renormalization used to 
treat the Coulomb interaction in the AGS equations is 
stretched to its limit of applicability in these regions of 
phase space characterized by small momentum transfer 
in the p- 12 G subsystem, a situation never encountered 
before in p-d [16] or d-a breakup [17]. 

For the n Be-proton test case, the picture is more com- 
plicated. For elastic scattering the two methods are in 
good agreement, however still fall short to describe the 
data. For the transfer cross section, CDCC constructed 
on the deuteron continuum underestimates the cross sec- 
tion compared to the solution of the AGS equations. And 
finally for the breakup observables, we only find good 
agreement between CDCC and Faddeev in certain re- 
gions of phase space, depending strongly on the choice 
of the basis used for the CDCC expansion. Specifically, 
in the energy regime dominated by p-n quasi free scat- 
tering (forward 10 Be angle) the representation based on 
the pn system (CDCC-TR*) accounts well for the full 
three-body effects, whereas for large angles of the de- 
tected fragment, corresponding to small excitations of 
the n— 10 Be system, the basis constructed from the n Be 
continuum (CDCC-BU) is more appropriate. A word of 
caution is required for the choice of the CDCC basis to be 
used for given kinematical regimes. It is also important 
to note that the rate of convergence of the CDCC observ- 
ables is very slow, particularly for CDCC-BU. This may 
in part explain the disagreement found at some angles 
and energies of the detected fragments. 

The CDCC equations attempt to produce a wavefunc- 
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tion that describes all the three-body effects, from small 
internal projectile distances to very large ones. In fact, 
all CDCC-BU observables are obtained here from the 
asymptotics of the three-body wavefunction. Since this 
is computationally very demanding, an alternative has 
been suggested that consists in using the wavefunction 
only in the range of the interactions, that is, inserted into 
the post form transition amplitude [42] . Further work on 
this topic is needed to explore such possibility. 



may be useful in practical calculations. It was shown in 
Refs. [16, 28, 43] that the R — > oo limit in Eq.(6b) exists 
for both parts separately and that the longer-range part 
of the breakup amplitude is given by 



B£\Z) = [1 + T pR (Z)G (Z)]v a [l + GW(Z)TZ%(Z)] 

(A2) 
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where p is the neutral particle and T pR {Z) — w pR + 
w pR Go{Z)T pR (Z) is the two charged particle screened 
Coulomb transition matrix. This part of the ampli- 
tude, called pure Coulomb breakup term, requires larger 
screening radius for convergence if Afc is small, but is 
simpler to calculate than the full U^ R \z). In contrast, 
the convergence with R is faster for the shorter-range 
part [utf? (Z) - B^ R) (Z)], even when Afc is small. There- 
fore in kinematical configurations of d+ 12 C breakup with 
small Afc it is sufficient to calculate [U^(Z) - B { Q R) (Z)} 
with standard parameters described in Sec. IV A, but the 
remaining term £?Q a (Z) needs considerably larger R. In 
the latter case we use the form 



APPENDIX A 

The kinematical situations characterized by small mo- 
mentum transfer Afc in the subsystem of charged parti- 
cles are sensitive to the screened Coulomb potential at 
large distances. In elastic scattering Afc may even van- 
ish, but the problem is resolved by separating the long- 
range part of the amplitude (4) and explicitly performing 
the R — ► oo limit. In breakup the transition operator is 
a Coulomb-distorted short-range operator, Afc is always 
nonzero, and the R — > oo limit can be reached at finite 
R with sufficient accuracy. Nevertheless, the decompo- 
sition of the breakup operator (2b) into two parts with 
different range properties 

u£\z) = b£\z) + [u£\z) b£\z)] (ai) 



Bg>(Z) = w pR G (Z)T pR (Z) W^Gi R \Z)T^(Z) 
+ T pR (Z)G (Z)v a Gi R Hz)T^(Z) 

(A3) 



that is equivalent to (A2) on-shcll. The partial-wave con- 
vergence is slowest for the last term w pR — W^jj" which we 
therefore calculate without partial- wave expansion. Rea- 
sonably converged d+ 12 C breakup results in the present 
paper are obtained with the screening radius up to 60 fm 
for the pure Coulomb breakup term. Partial waves with 
p- 12 C orbital angular momentum I < 38 and total angu- 
lar momentum J < 100 arc included in the calculation of 
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